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ABSTRACT 

We study the properties of the diffuse light in galaxy clusters forming in a large hydrodynamical cosmological 
simulation of the ACDM cosmology. The simulation includes a model for radiative cooling, star formation 
in dense cold gas, and feedback by SN-II explosions. We select clusters having mass M > 10 14 A _1 M Q and 
study the spatial distribution of their star particles. While most stellar light is concentrated in gravitationally 
bound galaxies orbiting in the cluster potential, we find evidence for a substantial diffuse component, which 
may account for the extended halos of light observed around central cD galaxies. We find that more massive 
simulated clusters have a larger fraction of stars in the diffuse light than the less massive ones. The intracluster 
light is more centrally concentrated than the galaxy light, and the stars in the diffuse component are on average 
older than the stars in cluster galaxies, supporting the view that the diffuse light is not a random sampling of 
the stellar population in the cluster galaxies. We thus expect that at least ~ 10% of the stars in a cluster may be 
distributed as intracluster light, largely hidden thus far due to its very low surface brightness. 
Subject headings: cluster: evolution, X-ray masses, diffuse light 
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1. INTRODUCTION 

The presence of diffuse 'intracluster light' in galaxy groups 
and clusters is now well established; observations by several 
groups provide estimates of the fraction of diffuse light and its 
distribution, using different techniques (see Arnaboldi 2003 
for a review). The fraction of stars contained in this space- 
filling component seems to increase strongly with the den- 
sity of the environment: from loose groups (< 2%, Castro- 
Rodriguez et al. 2003; Durrell et al. 2003) to Virgo-like 
( 10%; Feldmeier et al. 2003; Arnaboldi et al. 2003) and rich 
clusters (~ 20% or higher; Gonzalez et al. 2000, Feldmeier et 
al. 2002; Gal- Yam et al. 2003). This correlation may repre- 
sent an important clue for understanding the mechanisms that 
produce intracluster (IC) light and drive its evolution in the 
cluster environment. 

Cosmological simulations of structure formation facilitate 
studies of the diffuse light and its expected properties. Du- 
binski (1998) constructed compound models of disk galaxies 
and placed them into a partially evolved simulation of cluster 
formation, allowing an evolutionary study of the dark mat- 
ter and stellar components independently. Using an empirical 
method to identify stellar tracer particles in high-resolution 
dark matter (DM) simulations, Napolitano et al. (2003) stud- 
ied a Virgo-like cluster, finding evidence of a young dynami- 
cal age of the intracluster component. The main limitations in 
these approaches is the restriction to collisionless dynamics. 

1 INAF, Osservatorio Astronomico di Pino Torinese, Strada Osservatorio 
20, 10025 Pino Torinese, Italy (murante@to. astro. it,arnaboldi@to. astro. it) 

2 Astronomisches Institut der Universitaet, CH-4102 Binningen, Switzer- 
land (ortwin.gerhard@unibas.ch) 

3 Dipartimento di Astronomia dell'Universita di Trieste, via Tiepolo 11, 
1-34131 Trieste, Italy (borgani, cheng@ts.astro.it) 

4 INFN - National Institute for Nuclear Physics, Trieste, Italy 

5 Dipartimento di Fisica Generale "Amedeo Avogadro", Universita degli 
Studi di Torino, Via Giuria 1, 1-10125, Torino, Italy (diaferio@ph.unito.it) 

6 Dipartimento di Astronomia, Universita di Padova, vicolo 
dell' Osservatorio 2, 1-35122 Padova, Italy (kdolag,tormen@pd.astro.it) 

7 Dipartimento di Astronomia, Universita di Bologna, via Ranzani 1, 
1-40127 Bologna, Italy (moscardini@bo.astro.it) 

8 INAF, Osservatorio Astronomico di Trieste, via Tiepolo 11, 1-34131 
Trieste, Italy (tozzi@ts.astro.it) 



In this Letter, we analyze for the first time the IC light 
formed in a cosmological hydrodynamical simulation includ- 
ing a self-consistent model for star formation. In this method, 
no assumptions about the structural properties of the forming 
galaxies need to be made, and the gradual formation process 
of the stars, as well as their subsequent dynamical evolution in 
the non-linearly evolving gravitational potential can be seen 
as a direct consequence of the ACDM initial conditions. It 
is therefore of immediate interest whether this theoretical for- 
mation scenario makes predictions for IC light consistent with 
observations. Using a large volume of 192 3 /T 3 Mpc 3 , we can 
furthermore study a statistically significant sample of clusters 
at z = 0, and analyze the correlations of properties of diffuse 
light with, e.g., cluster mass and X-ray temperatures. 

2. COSMOLOGICAL SIMULATIONS 

We analyze the large scale cosmological hydrodynamical 
simulation (LSCS) of a "concordance" ACDM model (O m = 
0.3, D A = 0.7, n h = 0.019 h~ 2 , h = 0.7, a & = 0.8) of Borgani et 
al. (2004, B04). It was carried out with the massively par- 
allel Tree+SPH code GADGET (Springel et al. 2001), us- 
ing 480 3 DM particles and as many gas particles. For the 
periodic cube of size 192/T 1 Mpc, the mass resolution was 
thus m dm = 4.62 x l0 9 hr 1 M Q and m gas = 6.93 x IO^Mq; 
the Plummer-equivalent softening length was e = 7.5 h kpc 
(at z = 0). Besides gravity and hydrodynamics, the simula- 
tion accounts for star formation using a sub-resolution multi- 
phase model for the interstellar medium (Springel & Hern- 
quist 2003), feedback from supernovae explosions (including 
the effect of galactic outflows), radiative cooling of the gas 
(assuming zero metallicity) and heating by a photo-ionizing, 
uniform, time-dependent UV background (Haardt & Madau 
1999). 

Clusters are identified at z = using a standard friends- 
of-friends algorithm, with a linking length of 0.15 times the 
mean DM inter-particle separation. We identified 117 clus- 
ters with M > 1O 14 /i _1 M0. Cluster centers were placed at the 
position of the most bound particle belonging to each group. 
X-ray temperatures and masses were evaluated at the radius 
R200, which encloses an average density of 200 times the crit- 
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ical density. A detailed study of the X-ray properties of our 
clusters, together with a full description of the numerical sim- 
ulation can be found in B04. 

This cosmological simulation shows an encouraging agree- 
ment with some of the most important observed X-ray cluster 
properties. However, there are also a number of discrepancies 
that remain unaccounted for, as in other comparable numeri- 
cal work. For instance, the observed radial temperature pro- 
files in the cluster centers are not reproduced, and the fraction 
of collapsed baryonic mass ("cold" gas and stars) appears still 
too large. While it is likely that the resolution of these prob- 
lems will require an improved treatment of the intracluster gas 
physics, nonetheless this simulation represents a useful tool to 
study the physical properties of the diffuse light. 

2.1. Presence of diffuse light 

For all clusters identified in the cosmological simulation, 
we analyze the stellar distribution in the cluster volume. To 
compare with surface brightness measurements, we compute 
the projected density of stars by integrating along a line-of- 
sight (LOS). We then extract the 2D-radial profile of each 
cluster centered on its most bound particle, computing density 
profiles in 100 shells for the star component from 0.05/? 2 oo to 
2/?2oo- A stacked profile is obtained by averaging the shell 
densities with the same radius R/R200 over all clusters. For 
each cluster, the stars in galaxies other than the cD or brightest 
cluster galaxies (BCGs) are not included. 

For the galaxy identification we use the publicly available 
package SKID (Stadel 2001 9 ) to identify self-bound gas and 
star particle groups within individual clusters. We selected a 
scale of 20/T 1 kpc, comparable to our physical force resolu- 
tion, as the typical SKID length-scale r. SKID groups to- 
gether particles lying near local maxima of the density field, 
as determined using DM, star and gas particles. Then those 
star and gas particles which have a total energy T + V > 0, 
where T is the kinetic energy and V the local gravitational en- 
ergy, are removed from each group. All particles in a sphere 
of radius 2r are considered for evaluating the gravitational en- 
ergy. We discarded groups having less than 32 star particles. 

We then look for a halo stellar component in the stacked 2D 
profile, following Schombert (1986), by checking whether the 
average radial surface brightness profile curves upward in a 
(p,R 1/a ) dia gram, over significant intervals in R x l a . We fit a 
Sersic law to the inner parts of the stacked surface brightness 
profile, in a range of radii from the center out to a radius where 
the surface density is about one third of the central value. The 
result is shown in Fig. ^ a light excess to the inner Sersic 
profile is evident for radii R/R200 > 0.18. The deviation at 
large R of the stacked surface brightness profile from the Ser- 
sic's law is interpreted as being due to an extensive luminous 
halo. In observed cDs this occurs at R ~ 50 -f- 80 kpc, which is 
smaller than our measured value. Previous studies of cluster 
evolution (Dubinski 1998, Napolitano et al. 2003) produced 
clusters with BCGs whose density profile followed a de Vau- 
coleurs' law at all radii. In our simulations, where gas, stars 
and DM particles are followed self-consistently during struc- 
ture formation, cD halos do form. 

According to Dressier (1979) and Kormendy (1980) the up- 
turn in the cD surface brightness (p,R x l a ) plot, indicating the 
additional luminous halo, would occur at the projected radius 
where the stars become unbound from the central elliptical 

9 http://www-hpcc.astro.washington.edu/tools/skid.html 




0.6 0.8 1.0 1.2 

(R/R 2 oo) (,/ " 



FIG. 1 . — Schombert-like analysis on the stacked 2D radial density profile 
(BCG + ICL) of clusters in the simulation (triangles). The light excess is 
evident at large cluster radii. We used an non-linear Levenberg-Marquardt 
fit to the function logE(r) = log£«. -3.33[(r/r c )'/" - 1] (solid line), with 
best-fit parameters logS,. = 20.80, r e = 0.005, a = 3.66. Also shown are the 
averaged 2D density profile of stars in galaxi es (do tted line) and in the field 
(dashed line), as determined by SKID (see If 12.11 . In the inserts, we show 
the same analysis for the most luminous clusters with T > 4 keV (left panel), 
and for less luminous ones with < T < 2 keV (right panel). The resulting 
best-fit parameters are respectively logS c = 16.47, r e = 0.11, a = 1.24 and 
logS e = 23.1 1, r e = 0.00076, a = 4.37. In the main plot and in the inserts the 
unit (S/i?2oo)'/ a refers to the a values given by each Sersic profile. 



galaxy and orbit in the cluster potential. The velocity disper- 
sion profile would then rise at the radius where the change of 
slope takes place. A similar effect is observed in NGC 1399, 
a nearby cD galaxy (Arnaboldi et al. 1994, Napolitano et al. 
2002). 

In a cosmological simulation, we have the phase space in- 
formation for all particles; thus, we can study the dynamical 
behavior of those particles which populate the outer halo of 
cD galaxies. We will refer to all star particles grouped to any 
substructures by SKID as "bound", while the others will be 
named "unbound". The distinction between stars bound to a 
galaxy and unbound stars is workable except for the cluster 
center, where the most bound IC stars and the stars of the cen- 
tral cD fall in the same part of phase-space. 

Once the stellar particles in a given cluster are flagged, we 
build up the global 3D radial profiles for the "bound" and the 
"unbound" stars. Then we group all the clusters in classes 
with different X-ray temperature, and derive the average 3D 
stellar density profiles for the "bound" and "unbound" com- 
ponents. Fig. [2] shows that the "unbound" stars have a shal- 
lower radial profile than the "bound" component, in the range 
of radii from 0.05/? 200 to 0.3/?2oo- These stars are responsible 
for the additional light detected in cD halos, and build up the 
diffuse light in our clusters. 

3. PROPERTIES OF DIFFUSE LIGHT 

The size of our simulated volume allows us to study the 
physical properties of the IC stellar light in clusters statisti- 
cally, and explore their dependence on the total mass and/or 
the X-ray temperatures of clusters. For all the clusters se- 
lected in our database, we evaluate the logarithmic slope of 
the 3D profiles in a range of radii around where R* is the 
cluster radius where p*. bound = P*. unbound- In the radial interval 
R*/l .5 <R< 1 .5 R*, the profile slopes for the diffuse stars are 
in the range (-1, -3), while the profile slopes for the "bound" 
stars are in the range (-3,-5). Here the diffuse component has 
a shallower slope than the "bound" stars belonging to the cen- 
tral cD. At larger radii, the "bound" stars are those in cluster 
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FIG. 2. — 3D radial density profiles of the bound (open triangles) and 
unbound (filled diamonds) stars, for clusters divided in X-ray temperature 
classes. From top to bottom, the lines are for: global average, T < 2 keV, 
2 < T < 3 keV, 3 < T < 4 keV, T > 4 keV. The last four couples of density 
profiles have been divided by 10 2 , 10 4 , 10 6 , 10 s for clarity. The radius where 
the bound component dominates shifts toward smaller fractions of the aver- 
age R200, while the temperature, and thus the mass of clusters, increases. We 
also show the 3D radial density profiles of one cluster compared with its re- 
simulation RIO (see text). Solid lines refer to RIO, dashed lines to the LSCS 
cluster. All density profiles for "unbound" and "bound" stars have been di- 
vided by 7 ■ 10' and 10 12 , respectively. 



galaxies, and they become dominant over the more centrally 
concentrated IC light, see Fig. [1] and Fig. 

We also investigated whether any trends are present with 
cluster mass, either for the slope or for but no significant 
trend was detected. 

Another relevant quantity is the fraction / of diffuse vs. to- 
tal mass in stars: this quantity is a function of the cluster mass, 
with the more massive clusters having a larger fraction of IC 
light, see Fig. [3] We found that the amplitude of the corre- 
lation between the fraction of diffuse vs. total stellar mass 
depends on the value of the SKID length parameter r, but the 
slope of this correlation is almost unchanged when r changes. 
This is connected with the difficulty of separating the stars in 
the cD from the IC stars near the cluster center. Large values 
of r increase the number of stars in galaxies, but the fraction 
of unbound stars never drops to zero. In the most extreme 
case we tested (r = 40/T 1 kpc), only few of our 1 17 clusters 
drop to f < 10%. These objects are in the low-mass range 
(M<2-10 14 /z- 1 M Q ). 

We checked that r = 20/r'kpc ensures a clear dynami- 
cal separation between "bound" and "unbound" stars. The 
"unbound" stars have a 3D velocity dispersion comparable 
to that of the DM particles, while "bound" particles have a 
smaller velocity dispersion. Typical values for the most mas- 
sive cluster in our LSCS are er U nbound — &dm ~ 3000 kms" 1 , 
Abound » 800 kms -1 at R = 100/T 1 kpc. 

The simulation also records the age of formation of each 
star particle. In Fig. [5] we investigate the age distribution 
of the IC stellar components: stars in the diffuse component 
formed from their parent gas particles at an earlier average 
redshift (z rs 1.9), than "bound" stars (z w 1.7), with no ev- 
idence for a dependence on cluster mass. Our cosmologi- 
cal simulation predicts that the stars in the IC component are 
older. 

Two numerical effects can influence our results: the numer- 
ical resolution, and the parameter r for substructure identifi- 
cation. To study the first effect, we simulated a cluster with 
total mass 2.9 • 10 14 Ii~ 1 Mq, and amass resolution increased by 
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FIG. 3. — Left: Fraction of stellar mass in diffuse light vs. cluster mass. 
Dots are for clusters in the simulated volume; asterisks show the average 
values of this fraction in 9 mass bins with errorbars. Right panel: histograms 
of clusters over mean formation redshift, of their respective bound (dashed) 
and IC star particles (solid line). Mean formation redshifts are evaluated for 
each cluster as the average on the formation redshift of each star particle. 



a factor 3 and 10 respectively (R3 and R10 runs, here below). 
Softening was rescaled as m" 1 / 3 . Initial conditions for these 
simulations were generated using the "zoomed initial condi- 
tions" technique (Tormen et al. 1997), which increases the 
resolution in the cluster Lagrangian region while maintaining 
a coarse sampling of the surrounding structures to account for 
their tidal field. We find that the amplitude and slopes of the 
density profiles for both bound and unbound stars are almost 
unchanged at radii as large as w 200/;"' kpc (see Fig. [2). The 
lower resolution in our simulation may lead to an enhance- 
ment of the unbound population because of numerical over- 
merging. However, when we increase the resolution by a fac- 
tor 10 in mass and m 3 in force, the fraction / does not change 
substantially: we find / = 0.41 in LSCS, / = 0.43 for R3, and 
/ = 0.38 for R10 for the same cluster, suggesting that the nu- 
merical resolution has only a small effect. Moreover, a recent 
analysis of simulations of clusters having similar resolution 
to our R3 (Sommer-Larsen et al. 2004), has independently 
confirmed the / value reported here. 

The SKID length parameter does influence the behaviour of 
both "bound" and "unbound" components in a number of ex- 
pected trends, i.e., the mass of the cD galaxy and the value of 
increase. However, even when r = 40/T 1 kpc, the fraction 
of IC light is not smaller than 0.1 in the less massive clus- 
ters. Among those r values that we checked, r = 20/T 1 kpc 
ensures the best dynamical separation between the two stellar 
components at the cluster centers. 

4. DISCUSSION 

We used a cosmological simulation of (192/i _1 Mpc) 3 to 
study the statistical properties of the IC light in clusters of 
galaxies and the dependence of its physical properties on clus- 
ter mass and X-ray temperature. These predictions can be 
tested against known properties of cD halos and used to plan 
observational tests to understand the physical properties of IC 
light. 

The presence of the IC component is evident when the 
whole distribution of stars in the simulated clusters is anal- 
ysed in a way similar to Schombert's (1986) photometry of 
BCGs. Galaxies at the center of our simulated clusters have 
surface-brightness profiles which turn strongly upward in a 
(p,R l / a ) plot. This light excess can be explained as IC stars 



4 



Murante, G. et al. 



orbiting in the cluster potential. Integrating its density distri- 
bution along the LOS, the slopes from our simulations are in 
agreement with those observed for the surface brightness pro- 
files of the diffuse light in nearby clusters. In the Coma clus- 
ter, Bernstein et al. (1995) parametrize the surface brightness 
as r 13 and find that the diffuse light is best fit by /? = — 1.3 ±0.1. 
In the Fornax cluster, the surface brightness profile of the cD 
envelope of NGC 1399 follows a power law of the form oc r 13 
with /3 = -1.5 (Bicknell et al. 1989). 

At large cluster radii, the surface brightness profile of the 
IC light appears more centrally concentrated than the surface 
brightness profile of cluster galaxies (see Fig.[Qand Fig.|2j. 
From the simulations we also obtained the redshifts z/orm at 
which the stars formed: those in the IC component have a 
Z/orm distribution which differs from that in cluster galaxies, 
see Fig. [3] The "unbound" stars are formed earlier than the 
stars in galaxies. The prediction for an old stars' age in the 
diffuse component agrees with the HST observation of the 
IRGB stars in the Virgo IC field, e.g. t > 2Gyr (Durrell et 
al. 2002), and points toward the early tidal interactions as the 
preferred formation process for the IC light. The different age 
and spatial distribution of the stars in the diffuse component 
indicate that it is a stellar population that is not a random sam- 
pling of the stellar populations in cluster galaxies. 

The more massive clusters have the largest fraction / of 
diffuse light (Fig. [3}- It is /> 0.1 for cluster masses M > 
1O 14 /i _1 M0. Our simulations may thus explain the low in- 
ferred star-formation efficiency in clusters vs. less massive 
structures (David 1997). If only the bound stellar mass is ac- 
counted for in the ratio of the total cluster stellar mass vs. 



cluster gas mass in our LSCS, then this ratio decreases from 
groups to rich clusters. The observational trend would then 
be reproduced in the simulation. Similarly, the disagreement 
found between the amount of stars produced in clusters in our 
LSCS and in observed clusters (see B04) is less severe, if an 
IC component is present in real clusters and has been system- 
atically neglected when evaluating their internal stellar mass 
budget. 

The main result of this work is that large cosmological hy- 
drodynamical simulations are in qualitative agreement with 
the observed properties of diffuse light in galaxy clusters. A 
quantitative assessment will require additional numerical ef- 
forts and more observations. A detailed study of the dynami- 
cal history of the unbound stellar population in our simulation 
will be presented in a forthcoming work. 
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